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Abstract 

We performed Monte Carlo simulations of systems of wedge-shaped objects formed from Gay- 
Berne ellipsoids joined to Lennard- Jones spheres. We studied two different wedge shapes, one 
more asymmetric than the other. The bend and splay flexoelectric coefficients were measured in 
the isotropic and smectic phases using linear response theory, and found to be negligibly small 
in the isotropic phase. We found a close connection between the properties of the intermolecular 
potential and the flexoelectric coefficients measured in the smectic phase. In particular, we found 
negligible bend coefficients for both shapes and a larger magnitude of the splay coefficient for 
the more prominent wedge, in accord with Meyer's original mechanism for flexoelectricity. The 
less prominent wedge produced a splay flexoelectric coefficient with the opposite sign due to the 
attractive tail of the intermolecular potential and the relative narrowness of the molecular head. 



1 Introduction 

In the flexoelectric effect [|l| a director deformation produces an electrical polarization, similar to 
the phenomenon of piezoelectricity in solid crystals. The flexoelectric effect was first proposed 
by Meyer ||l| who considered asymmetric molecules either wedge-shaped with longitudinal dipole 
moments or "banana" shaped with transverse dipoles. In the absence of a director deformation, the 
packing of the molecules is similar to that of ellipsoidal shaped or rodlike molecules (the additional 
asymmetries of the molecules have negligible effect) and the average polarization is zero. However, 
when a splay is imposed upon a system of wedges or a bend upon a system of bananas, the preferred 
packing of the molecules results in a net alignment of dipoles leading to an overall polarization of 
the medium (see figures ^ and ||). Alternatively, an applied electric field which aligns the dipoles 
induces a splay or bend in the appropriately shaped system — this is sometimes called the inverse 
flexoelectric effect. In either case, the net polarization P and the elastic deformations are related by 
the flexoelectric coefficients en and 633 introduced by Meyer through the following linear response 
relation: 

P = eiin(V-n) + e33nx (Vxn), (1) 

where fi is the director. The first term on the right hand side of this equation corresponds to 
splay flexoelectricity (relevant for wedges) and the second term to bend flexoelectricity (relevant 
for bananas). 

Subsequent to Meyer's work. Frost and Marcerou [^ proposed a flexoelectric mechanism based 
on molecular quadrupoles requiring neither the shape asymmetries nor the dipole moments of 
Meyer's original argument. Rather, molecular quadrupoles allow uneven charge distributions lead- 
ing to polarizations in any given volume when a splay is imposed (see figure ^). 

Both the Meyer's dipole and Frost and Marcerou's quadrupole mechanisms have been observed 
experimentally |^ Q, ^, |6|, ^. Typically, the flexoelectric coefficients are measured over a range of 
temperatures in the nematic phase and their variation compared with that of the nematic order 
parameter S. Although it has been shown that for both the dipole and quadrupole mechanisms 
there are contributions to the flexoelectric coefficients involving several powers of S, Marcerou and 
Frost used the dominant contributions — S for the quadrupole mechanism and 5^ for the dipole 
mechanism — in the analysis of their experimental data [Q. They found flexoelectric coefficients 
proportional to S for symmetric, non-polar molecules, clearly demonstrating the quadrupole mech- 
anism. Quadratic variation with S was seen for banana-shaped molecules with strong, transverse 
dipole moments. 

In the smectic A phase, an additional flexoelectric coefficient 622 arises, representing the coupling 
between the net polarization and variations in the smectic layer spacing g, g, ^ . Specifically, there 
is an additional contribution to the right-hand side of equation (||) of the form: 622 ^"t^, where 
u is the displacement of the smectic layers whose normals are parallel to the z axis. Frost and 
Fershan Q found this additional term to be negligible experimentally. The flexoelectricity of the 
smectic C phase involves a total of 14 coefficients and is discussed in |^]. 

A systematic method for calculating flexoelectric coefficients based on molecular shape and 
multipole properties would be very valuable. Mean-field theories pi, 11^, 11, 12, 13| do not consider 



the short-ranged fiuctuations in molecular alignment which can be quite important. Computer 
simulations offer a way to assess the molecular origins of flexoelectricity. In particular simulations 
can focus strictly on Meyer's packing ideas without the complication of dipolar interactions which 
could lead to antiparallel alignment of side-by-side molecules. To date only one simulation study 



TM of flexoelectricity has been carried out. In this study the flexoelectric coefficients were evaluated 



for wedge-shaped molecules interacting via a generalized Gay-Berne potential [15|. The molecules 



were modeled in this study (as well as in ours) by a Gay-Berne ellipsoid with a Lennard-Jones 



sphere attached near one end (see figure U). Two sets of parameters were considered in ref.|14], one 
with a slightly more pronounced wedge shape than the other. The flexoelectric coefficients were 
evaluated using microscopic expressions based on density functional theory, and a larger value of 
eii was obtained for the more pronounced wedge. The bend flexoelectric coefficient was nearly zero, 
in agreement with Meyer's suggestion that bend flexoelectricity should not appear in a system of 
wedges. 

In this paper we consider a similar model but extend the study in several ways. First, we 
evaluate the flexoelectric coefficients using linear response theory and the fluctuation-dissipation 
theorem, which provide a computationally simpler and more direct means of evaluation compared 
with density functional theory. Second, we consider two model sets of molecular parameters, one 
representing a molecule with significantly more asymmetry than either molecule considered in 



ref.|14]. By exploring the intermolecular potentials for our two sets of parameters we find that the 



less prominent wedges (similar to the less prominent ones considered in ref.|14]) prefer to align with 
their larger ends tilted toward each other, while the more prominent wedges prefer to align with 
their larger ends tilted away from each other. These opposing tendencies in turn lead to opposite 
signs of the splay flexoelectric coefficients, as one might expect from Meyer's model. As in ref. |14| 
we find a negligible value of the bend fiexoelectric coefficient. 

The outline of this paper is as follows. In the next section we present the details of our modeling 
of wedge-shaped molecules. Section III discusses the linear response theory used to measure the 
flexoelectric coefficients in our simulation. Our results are presented in section IV and we offer 
some concluding remarks in the flnal section. 

2 Molecular Shape Modeling 

Using an approach similar to that of [^J] we constructed a wedge-shaped molecule from a standard 
Gay-Berne ellipsoid (or rod) with a sphere added near one end (see figure^. The net interaction 
potential between two wedge-shaped molecules labeled 1 and 2 then consists of four terms, namely 

^tot — ^ rodl—rod2 ~r ^ spherel—sphere2 ~r l-'rodl—sphere2 ~r fJ spherel—rod2- V-^j 

where Urodi-rod2 is given by the original Gay-Berne potential [p^] : 

Urodl-rod2 (ui, U2, r) = 4e (Ui, U2, r) 

12 < 



X ,. ! ., - ,, ! ., , (3) 

Lr - a (ui, U2, r) -f- (To J Lr - cj (ui,U2,r) + CTo J J 

where ui, U2 give the orientations of the long axes of rods 1 and 2, respectively, and r = ri — r2, 
with the centers of the rods located at positions ri and r2. The parameter (T(ui,U2,f) is the 
separation between the rods at which the potential vanishes, and thus represents the shape of the 
rods. Its explicit form is 



0-(ui,U2,f) = (To 




r • ui -l-r • U2) 

1 + X(ui -112) 



(f • Ui - f • U2)^ 



1/2 

(4) 



1-X(ui -ua) 
where cTq = ax_ (defined below) and x is 

{(^||M)'-i}/{(^||/^±)' + i}- (5) 



X 



Here an is the separation between two rods when they are oriented end-to-end with Urodi-rod2 = 
0, and a± is the corresponding separation when the two rods are side-by-side. The well depth 
e (ui, U2, f ), representing the anisotropy of the attractive interactions, is written as 

e(ui,U2,r) = Soe" {ui,U2) e'f" {u-i_,U2,r) , (6) 

where 

e(ui,U2) = {l-x'(ui-U2f}''^', (7) 

and 



e'(ui,U2,f) = 1 - -X' 



1 , j (f -ui +f -112) 



2 



2^ [ 1 + X' (ill • ii2) 

1 - X (ui • U2) J 
with x' defined in terms of ey and e±, the end-to-end and side- by-side well depths, respectively, as 



V= l-CqiAx)"" / l+(qiAx)"'T (9) 



We measure all of our physical quantities in reduced units in terms of the energy scale Eq and the 
length scale (Tq. For the adjustable parameters appearing in equations (|6|), (^ and (^), we used the 



values originally proposed by Gay and Berne |1(:]: ^ = 2^u = 1, and £±/£\\ = 5. Our choices for 
the rod shape parameters (T\\ and a^ are discussed below. 

The interaction Uspherei-sphere2 is given by the ordinary Lennard-Jones potential: 



dV^ fd 



Uspherel-sphere2 (r) = 4eo \ [~) ~ ( ~ ) f ' (^0) 



where d is the separation between the two spheres at which the potential Uspherei-sphere2 vanishes; 
i.e., it is a measure of the diameter of the sphere. The relative position vector r is measured from 
the center of sphere 1 to the center of the sphere 2. 

The interaction between the rodlike part of one molecule and the sphere of the other molecule, 
Urodi-sphere2i is given by a Gay-Berne potential generalized to mimic the interaction between 



nonequivalent particles [15|. The range parameter ars (ui, Y12) and the energy parameter e^s (ui, Uj) 
used in this potential are generalizations of the corresponding parameters in equations (H) and (pi), 
namely [|l5| 

a,.(ui,f) = a-(l-Xr.(r-Cii))-i/2 (11) 



and 

£rs{<il,r) = 1, (12) 

with Xrs defined as 



a2-a2 



and 

a- = ^{al + d')'/'. (14) 

In the rod-sphere potential the relative position vector r is measured from the center of the rod to 
the center of the sphere. 

We performed simulations for two sets of molecular parameters: a slightly wedge-shaped object 
with a± = 1.0, an = 2.6, d = 0.94, and a more prominent wedge with parameters a± = 1.0, fJii = 
2.4, d = 1.3. In the former case the center of the sphere was located at a distance D = 1.3 from 
the center of the rod, and at a distance D = 1.2 in the latter case. In terms of the steric dipole 
moment p* = (47r/3)(d/2)^Z) introduced in [14|, our parameters correspond to values of p* of 0.565 
and 1.38 respectively (in [|l4| systems with dipole moments of 0.524 and 0.662 were studied). 

To ascertain the effective shape of the composite rod-sphere molecule we computed equipotential 
contours (with contour values close to zero). The contour plots are shown for the two sets of 
parameters in figures ^ and ^ clearly indicating a wedge-like shape. To explore the interaction and 
local packing of the molecules we computed the depths of the potential wells of Utot for several 
different relative orientations of a pair of molecules with relative tilt 9. The results for the two sets 
of molecular parameters are shown in figures and|8|. For the purposes of comparison corresponding 
curves for the original Gay-Berne potential Urodi-rod2 are shown in figure ^ Comparing figures ^ 
and g with figure P, we note that in the splay configuration the wedges do not exhibit the left-right 
tilt symmetry (9^6 + 180°) seen in the Gay-Berne case. This result is not surprising: tilting the 
larger end of one wedge away from the larger end of another wedge (corresponding to ^ < 180° in 
the figures) should yield a different energy than tilting the larger ends toward each other. More 
importantly, for each set of wedge parameters there is an absolute minimum in the potential energy 
corresponding to a finite but small splay angle. For the more prominent wedge shape (figure |8|) 
this angle is approximately 10°, i.e., the pair of wedges prefer to align with their larger ends tilted 
away from each other. For the less prominent wedge shape (figure 0) the angle corresponding to 
the absolute minimum is nearly 360°, i.e., the wedges prefer to align with their larger ends tilted 
toward each other. As we shall see below this difference yields flexoelectric coefficients of opposite 
signs for the two molecular shapes. While it might seem surprising at first glance that the less 
prominent wedges tilt toward each other, given the attractive interaction between two spheres and 
the fact that d < a± in this case it is a reasonable result. For the more prominent wedges the larger 
repulsive core of the spheres leads to the expected tilting of the larger ends away from each other. 

While the absolute minima of the pair potentials for the two sets of wedge parameters correspond 
to nonzero splay, note that this state is only slightly preferred over the aligned state, 9 = 0°, for 
the parameters of figure and the antiparallel aligned state, 9 = 180°, for the more prominent 
wedge of figure ^. There is also a substantial subsidiary potential minimum for antiparallel slightly 
wedge-shaped molecules. In our simulations of a system of 256 wedge-shaped molecules using both 
sets of molecular parameters we found no net spontaneous splay or electric polarization (assuming 
molecular dipole moments parallel to the long axis of the wedge). However, if an electric field 



were applied to the system of prominent wedges, we would expect the potential well corresponding 
to nonzero splay to become deeper relative to that of the antiparallel state 9 = 180°, leading to 
splay flexoelectricity. For the slightly wedge-shaped objects of figure |^ the shallowness of the well 
corresponding to nonzero splay makes it less obvious that splay flexoelectricity will exist. However, 
in our simulations we do in fact find splay flexoelectricity in this case, albeit smaller than in the 
case of the more prominent wedges and with the opposite sign. 

We now turn to a discussion of how we extract the flexoelectric coefficients from our simulations. 

3 Calculating flexoelectric coefficients 

To measure flexoelectric coefficients in our simulation, we used the linear response theory of Nemtsov 



and Osipov |17]. The flexoelectric coefficients in this method are related to the response function 
of the system to an orientational stress. Using the fluctuation-dissipation theorem the response 
function can be found from correlation functions of the polarization density and the orientational 
stress tensor. The latter tensor is conjugate to the orientational strain which yields flexoelectricity. 
Thus, a calculation of the relevant correlation functions yields the flexoelectric coefficients. 

Specifically, the splay and bend flexoelectric coefficients in the Nemstov-Osipov formalism are 
given by 

eii = — £"0/3^6^/3^ n^n^/ 2 

633 = Eai3-yeai3,in^nn/2, (15) 

where we use the summation convention over the Greek indices (summed over the coordinate direc- 
tions X, y and z). The tensor Bq,^^ is the antisymmetric Levi-Civita tensor, while the antisymmetric 
tensor Eap^j is the response function satisfying 

Pa = -E'a/377/37 (16) 

where ■jais is the orientational strain tensor given by 

dxp 

Here 9a denotes the rotation angle of the director about the coordinate axis labeled by a. For small 
director deformations 5h, 9 ~ sin ~ n x 5h, which then yields 

7/37 = efSf^un^,—. (18) 

Symmetry considerations show that Ea/3-y has four independent components in the nematic and 
smectic A phases. 

Using the fluctuation-dissipation theorem the components of the response function E^fB-y are 
given by correlation functions of the orientational stress tensor and polarization: 

Eaf3^ = -y (T^fS^'Pa) ■ (19) 



Here (3 = l/ksT, V is the volume of the system, V = J2i l^ii it^i is the dipole moment of molecule 
i), and tTq^ is the static orientational stress tensor given by 

^a/3 = 2 ^ '^ijpTija- (20) 

In the equation above Vij is the relative position vector of molecules i and j, and tjj is the torque 
exerted by molecule j on i. 

In our simulations, then, we calculate the components E^is^ of the response function by com- 
puting the correlation functions in equation (jl^) as time averages (or MC cycle averages) and then 
calculate en and 633 from equation ( [l5| ) (the director is also computed during the simulations). 
We set the magnitude of the molecular dipole moment to unity; its direction is given by u, the long 
axis of the wedge. The torque on molecule i due to our generalized potential Utot is given by ||l| 



Ti = Yl '^ij = ^i ^ i-^uUtot) ■ (21) 

4 Results 

We performed constant temperature and pressure (with P* = Pa^/so = 10) Monte Carlo simula- 
tions on systems of 256 wedges using the two sets of molecular parameters described in Section II. 
The systems were equilibrated for 250,000-500,000 cycles and cooled in dimensionless temperature 



steps of 0.1 (the dimensionless temperature is defined by T* = kBT/so)- As in [14| we found a 
strong preference for the system to form a smectic A phase, even when we removed the attractive 
part of the potential (in the ordinary GB system, removal of the attractive forces tends to stabilize 
the nematic phase [jl9|). The use of a different set of Gay-Berne parameters for the attractive 
portion of the potential [£0| (which enhance the stability of the nematic phase for a system of 
GB ellipsoids) did not stabilize the nematic phase in the present case of wedge-like molecules — 
another indication that the repulsive core shape is the dominant factor in the possible formation 
of a nematic phase. 

Stelzer et al. [14| found that the nematic phase is absent for dipole moment p* = 0.814, though 
they were able to produce nematic phases over narrow temperature ranges for p* = 0.524 and 
0.662. Thus, our inability to produce a stable nematic phase for our molecules with p* = 1.38 is 
consistent with [|lj], but our inability to produce a stable nematic for our less asymmetric shape 
with p* = 0.565 is not, and the reason for this discrepancy is not clear. We note that experimentally 
the splay and bend flexoelectric coefficients in the nematic and smectic phases have been found 
to be virtually identical Q, while the fiexoelectric coefficient 622 which appears in the smectic A 
phase but not the nematic is negligibly small. Thus, we proceeded to measure the splay and bend 
flexoelectric coefficients in the smectic phase using the method outlined in section III. 

Interestingly, the smectic layers showed a distinct domain structure (see figure p!oD characterized 
by alternating regions of parallel molecular alignment (compare with the Gay-Berne smectic also 
shown in figure |l^). This structure seems consistent with the idea that the parallel wedges prefer 
splay; splay occurs in each domain while alternation of the alignment between domains maintains 
the fiat, parallel smectic layers. 

Data for the small (p* = 0.565) and large {p* = 1.38) wedges and for Gay-Berne ellipsoids 
are shown in Table Q. The Gay-Berne ellipsoids interact via the rod potential equation (^) with 



parameters a± = 1.0, an = 3.0,^ = 2, and i/ = 1. We note that the values of the flexoelectric 
coefficients for the system of ellipsoids are zero to within the computed error; thus in accord with 
Meyer's ideas this system of symmetric objects does not exhibit flexoelectricity. In the isotropic 
phase both the large and small wedges also exhibit no flexoelectricity. However, in the smectic 
phases of both wedges splay flexoelectricity appears. Note the difference in sign between the en 
data in the smectic phase for the large and small wedges, consistent with their opposite splays. 
The actual values for the signs also are as expected: considering again figure |l] with the director fi 
taken to point upwards, the splay shown is then positive and the resulting polarization is parallel 
to n, implying that en > (recall equation (||)); this is the case for the large wedge. For the 
small wedge, the splay is negative but the polarization is still parallel to n, implying en < as 
observed. The magnitude of en is also larger for the large wedge consistent with Meyer's excluded 
volume mechanism for flexoelectricity. The average values for ess, the bend flexoelectric coefficient, 
are much smaller in the smectic phases than the corresponding en values and are zero to within 
the computed error. We also note that the bend coefficients varied in sign over the course of the 
simulation whereas the values of en maintained a consistent sign during the runs. 

5 Conclusion 

By simulating a system of wedge-like molecules formed from Gay-Berne ellipsoids and Lennard- 
Jones spheres we have explored some of the molecular origins of flexoelectricity. We measured both 
the bend and splay flexoelectric coefficients using linear response theory, which yields the coefficients 
in terms of correlation functions of the molecular torque and orientation vector. We studied wedges 
with two different parameterizations, and found a close connection between the properties of the 
intermolecular potential and the flexoelectric response of the system. In particular, wedge-shaped 
molecules do not produce a measurable bend flexoelectric coefficient, and a more prominent wedge- 
shaped object produces a larger splay flexoelectric coefficient, in accord with Meyer's original ideas 
on the origins of flexoelectricity. In the case of the less prominently shaped wedge we have obtained 
a splay flexoelectric coefficient with the opposite sign to that of the more prominent wedge due 
to the attractive tail of the intermolecular potential and the relative narrowness of the molecular 
head. 
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Figure 1: (a) Wedges with longitudinal dipoles under normal nematic conditions. There is no splay 
and no net polarization, (b) Under an applied splay, the preferred wedge alignment results in a 
net polarization. Alternatively, an applied field induces a splay due to the wedge shape of the 
molecules. 

Figure 2: (a) "Bananas" with transverse dipoles under normal nematic conditions. There is no 
bend and no net polarization, (b) Under an applied bend, the preferred banana alignment results 
in a net polarization. Alternatively, an applied field induces a bend due to the banana shape of the 
molecules. 

Figure 3: Quadrupoles with a splay imposed. Within the central "layer", plus charges from above 
are allowed to enter, while plus charges below are expelled, leading to a net polarization upwards. 

Figure 4: Illustration of the basic geometric parameters of the wedge-shaped molecule composed 
of a rod and sphere. 

Figure 5: Potential energy contours showing the shape of the rod-sphere composite with parameters 
a± = 1.0, cr|| = 2.6, d = 0.94: (a) side view showing wedge-like asymmetry and (b) top view 
showing axial symmetry. The top view contours are not completely circular due to a finite number 
of sample points. The side view data was generated by fixing one wedge in the yz plane as shown, 
while a second wedge pointing up out of the plane and with its narrow end just touching the plane 
"scanned" across the first wedge. For the top view data, both wedges were parallel; one remained 
fixed while the other was moved around the first in the same plane. 

Figure 6: Potential energy contours showing the shape of the rod-sphere composite for a more 
prominent wedge with parameters a± = 1.0, a\\ = 2.4, d = 1.3. 



Figure 7: Minimum well depths (calculated as a function of molecular separation) for molecular 
parameters a± = 1.0, crii = 2.6, d = 0.94 for various relative orientations as shown: (a) splay, (b) 
twist and (c) bend configurations. 

Figure 8: Same as figure |^ but for the more prominent wedges shown above (figure |6|) with param- 
eters a± = 1.0, o-|| = 2.4, d = 1.3. 

Figure 9: Same as figure but for Gay-Berne molecules with parameters (t± = 1.0, au = 3.0. 

Figure 10: View of the molecular configurations in the smectic layers for (a) wedges and (b) standard 
Gay-Berne ellipsoids. The arrows indicate the direction of the orientation vectors u whose heads 
correspond to the wide ends of the wedges. The domains in the wedge case can be clearly seen. 
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Table 1: Flexoelectric coefficients for the small {p* = 0.565) and large (p* = 1.38) wedges and for 
Gay-Berne ellipsoids. Data for the small wedge was obtained in the isotropic phase at T* = 3.5 
and in the smectic phase at T* = 2.9; corresponding values for the large wedge were 2.5 and 1.9. 
Data for the Gay-Berne ellipsoids was obtained in the smectic phase at T* = 0.745 and P* = 2.5. 
The nematic order parameter in each of the smectic phases was 0.9. 





en 


633 


small wedge, isotropic phase 

large wedge, isotropic phase 

Gay-Berne ellipsoids, nematic phase 


0.367 ±0.522 

-1.39 ±1.15 

-0.394 ± 0.581 


-0.394 ± 0.383 

0.148 ±0.251 

-0.0835 ±0.117 


small wedge, smectic phase 
large wedge, smectic phase 


-2.08 ±0.211 
13.6 ± 0.052 


-0.061 ±0.054 
-0.005 ±0.012 
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